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Abstract 

We consider a relativistic, degenerate, electron gas under the influence of a 
strong magnetic field, which describes magnetized white dwarfs. Landau quanti- 
zation changes the density of states available to the electrons, thus modifying the 
underlying equation of state. In the presence of very strong magnetic fields a max- 
imum of either one, two or three Landau level(s) is/are occupied. We obtain the 
mass-radius relations for such white dwarfs and their detailed investigation leads us 
to propose the existence of white dwarfs having a mass ~ 2.3-M©, which overwhelm- 
ingly exceeds the Chandrasekhar mass limit. 
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Introduction 

Chandrasekhar first showed the maximum possible mass for a white dwarf (WD) to be 
~ 1.44M pp. Later, several magnetized WDs have been discovered with surface fields 
of 10 5 — 10 9 G [2], [3], [5], [6]. Anticipating that the central field, which cannot be 
probed directly, will be much stronger than that on the surface, Ostriker and Hartwick 
pp constructed models of WDs with B ~ 10 12 G at the center but with a much smaller 
field at the surface. Recent observations of peculiar Type la supernovae - SN2006gz, 
SN2007if, SN2009dc - seem to suggest super- Chandrasekhar mass WDs as their most 
likely progenitors [8]. These WDs, found in binary systems, stabilize by rotating and 
accreting matter [5J. In this essay we propose a mechanism by which the WDs exceed 
the Chandrasekhar mass limit in presence of strong magnetic fields at their centers. In a 
simple theoretical framework existence of such stars has been reported recently [10J. 

Here we consider strongly magnetized WDs having interior magnetic field > 10 15 G. 
These Landau quantized electronic systems have at the most one, two or three occupied 
Landau level(s) (LL(s)). We obtain an exciting possibility of a mass of WDs ~ 2.3M . 



Equation of state and solution procedure 

Basic equations 

The Landau quantized energy states of a free electron in a uniform, static, magnetic field 
are given by [UJ 



where v denotes the LL, given by 



v = j + 2 + °, ( 2 ) 



when j being the principal quantum number, a — ±|, p z the momentum of the electron 
along the z-axis, B the magnetic field, B c a critical magnetic field defined by 

B c = m 2 e c 3 /he = 4.414 x 10 13 G, (3) 

where e is the charge of the electron, c the speed of light, m e the rest-mass of the electron, 
h the Planck's constant. Our interest is to study the effect of B > B c on the relativistic, 
degenerate electron gas. 

The Fermi energy of electrons in units of m e c 2 for a given v is given by 

el = XF (v) 2 + (l + 2v^ r ), (4) 

where Xp{y) is the Fermi momentum in units m e c. 
As xp{v) 2 > 0, the maximum number of occupied LLs 

\ D / nearest lowest integer 



2 



where Bp = B/B c , ep max the dimensionless maximum Fermi energy of the system. 
Following [TT] we write the electron number density 



v / e v= q 

where the Compton wavelength of the electron A e = h/m e c; the matter density 

p = fi e m H n e , (7) 

where /i e is the mean molecular weight per electron (which we choose to be 2) and thh 
the mass of hydrogen atom; the electron degeneracy pressure 



where 

rj(y) = \y^i + y 2 - \ Hv + V^ + y 2 )- (9) 



Procedure 



From equation (jSJ) we see that if Bd 3> 1, then v m is small, implying electrons are 
restricted to the lower LLs. Since we are interested in this regime, depending on the 
specific values of B D and eF max , we have one- level (0 < v m < 1), two-level (1 < v m < 2) 
and three- level (2 < v m < 3) systems. By eliminating Xp(i/) numerically from equations 
(J7J) and (jS|), we obtain the P — p relation which is the equation of state (EOS). Figure 1 
shows EOSs for the different cases given in Table 1. 



Table 1 

Parameters for the equations of state in Figure 1. 



^Fmax 




B D 


B in units of 10 15 G 
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0.033 
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0.5 


0.022 




1 


199.5 


8.81 


20 


2 


99.75 


4.40 




3 


66.5 


2.94 




1 


19999.5 


882.78 


200 


2 


9999.75 


441.38 




3 


6666.5 


294.26 



The hydrostatic equilibrium of the WD in presence of constant magnetic field is de- 
scribed by [12] 

1 d fr 2 dP\ 



r 2 dr \ p dr ) ^' ^ ^ 



with the boundary conditions: 



p(r = 0)= Pc (11) 

and 

where p c is the central density of the WD. 

The radius R of the WD is obtained by solving equation ( FlOi) and is the value of 
r where density goes to zero. The mass M of the WD, which is approximated to be 
spherical, is obtained by integrating the following equation: 

dM „ 

__ = 4vrr 2 p. 13 
dr 



Constructing white dwarfs 

Equations of State 

Let us consider the case ep max = 20. The solid curve in Figure 1(b) is free of any kink 
(only ground LL occupied). The dotted and dashed curves have one kink (ground and 
first LLs occupied) and two kinks (ground, first and second LLs occupied) respectively. 
Any kink represents a transition from a lower to upper LL. 

Mass-radius relations 



Table 2 

Parameters for the fitting function for the equation of state shown in Figure 1(d). 



Pd in units of 2 x 10 9 gm/cc 


r 


K in CGS units 


- 0.096 


2.9 


4.055 


0.096 - 0.307 


2.4 


1.284 


0.307 - 1.128 


2.1 


0.914 


1.128 - 2.117 


0.35 


1.105 


2.117 - 2.956 


4/3 


0.522 


2.956 - 3.842 


2.0 


0.246 


3.842 - 4.651 


0.35 


2.225 


4.651 - 6.116 


4/3 


0.496 


6.116 - 7.37 


2.0 


0.142 



Figure 2 shows the mass-radius relations, where each point in the curves corresponds to 
a WD with a particular value of p c which is supplied as a boundary condition. 

In Figure 2(a) we see that as p c increases, both Mp and Ru increase and then at higher 
values of p c , Rd becomes nearly independent of Md- We note that the most massive WD 
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Figure 1: Equations of state in a strong magnetic field (given in Table 1) for (a) ep Ti 



(b) t Fri 



20, (c) e Fr 



200. In all three cases the solid, dotted and dashed lines 



indicate one-level, two-level and three-level systems respectively. In (d) the solid line is 
same as the dashed line in (b), but fitted with the dotted line by analytical formalism 
(see text for details). Here Pd is the pressure in units of 2.668 x 10 27 erg/cc and pn is 
the density in units of 2 x 10 9 gm/cc. 



on this curve has a mass ~ 2.3 M , which corresponds to the maximum density point 
of the solid curve in Figure 1(b). This denotes the density at which the ground LL is 
completely filled. Thus a WD with this p c and a magnetic field strength of B = 8.81 x 10 15 
G has a mass greater than the Chandrasekhar limit (for details see [13]). 

The turning point in Figure 2(b), after attending the maximum mass ~ 2.3M , corre- 
sponds to the kink in the corresponding EOS given in Figure 1(b). During the transition 
from ground to first LL, the radius and mass both decrease with increasing p c . Then 
there is a brief range of densities where the mass increases as the radius decreases and 
ultimately at very high densities the radius is nearly independent of the mass. 

In Figure 2(c) we see a repetition of the features in Figure 2(b) twice, which corre- 
sponds to the two kinks in the corresponding EOS in Figure 1(b). (The maximum mass 
~ 2.3M Q occurs at that p c where the ground LL is completely filled and transition to the 
first LL is about to start.) 
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Figure 2: Mass-radius relations with ep-max = 20 for (a) one-level system, (b) two- level 
system, (c) three-level system. Here Md is the mass of the white dwarf in units of M & 
and Rd is its radius in units of 10 8 cm (the solid, dotted and dashed lines have the same 
meaning as in Figure 1). 



For comparison let us recall the Lane-Emden solution for the classical non-magnetic 
case, which assumes a polytropic EOS 

P = Kp r , (14) 
giving rise to the following relations [12]: 

R oc pT 1 (15) 

and 



M oc p c 2 . (16) 

This means that if T > 2, R increases with p c , if F — 2, R is independent of p c , if T > 4/3, 
M increases with p c and if T = 4/3, M is independent of p c . This is exactly what is 
observed in Figure 2, if the EOSs given in Figure 1 are fitted by the polytropic EOS, 
adopting constant values of T in different density ranges (see Table 2). For instance in 
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Figure 1(b) (dashed line) and Figure 2(c), at very low densities, T = 3 and up to the 
turning point density (the kink in the EOS) the radius keeps increasing with mass and 
then becomes nearly independent of mass when T ~ 2. Then T suddenly drops to a small 
value ~ 0.35 which marks the onset of the transition from ground LL to first LL. In this 
region the pressure becomes independent of density, revealing an unstable zone in the 
EOS. As the density increases further, T approaches the relativistic value of 4/3. In this 
regime we see that the radius decreases slightly as the mass does not change significantly. 
Next T takes up a value of 2 and the above phenomena will be repeated. 



(a) 



(b) 





Pd ivi d 

Figure 3: Comparison with Chandrasekhar's non-magnetic results for ep m ax — 20. (a) 
Equations of state - the solid line represents Chandrasekhar's equation of state. The dot- 
dashed, dotted and dashed lines represent the one-level, two-level and three-level systems 
respectively. The equation of state for u m = 20 is also shown, which appears as a series 
of kinks on top of the solid line, (b) Mass-radius relations - the vertical line marks the 
1.44M limit and the solid line represents Chandrasekhar's mass-radius relation. From 
top to bottom the other lines represent u m = 500 , 20 , 3 , 2 and 1 respectively (the y-axis 
is in log scale). 



Discussions 

In Figure 3 we put together our results along with that of Chandrasekhar. Figure 
3(a) clearly indicates that as the magnetic field decreases, the EOS approaches Chan- 
drasekhar's EOS. 

Figure 3(b) represents the corresponding mass-radius relations. We observe that as 
the magnetic field decreases, the mass-radius relation approaches that obtained by Chan- 
drasekhar. We also note that as the magnetic field increases the WDs become more 
and more compact in size and the probability that they will have masses exceeding the 
Chandrasekhar limit also increases. 

Figure 4 shows the variation of mass as a function of density within a magnetized WD 
for three different magnetic fields. In all the cases, we note that, by the time the density 
falls to about half the value of p c , the mass has already crossed the Chandrasekhar limit 
as indicated by the horizontal line. Hence the effect of magnetic field is restricted to the 
high density regime, where the field remains essentially constant even in reality. This 
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Figure 4: Mass as a function of density inside a magnetized white dwarf with Bo (a) 
199.5, (b) 99.75, (c) 66.5, with ep max = 20. The horizontal line indicates the 1.44M 
limit. See Table 1 for details. 



justifies our choice of constant field. 



Conclusions 



We summarize the important findings of this work as follows: 

• A transition from the lower to upper LL represents a kink in the EOS. If p c lies at 
the kink appearing at the ground to first LL transition, the WD has the maximum 
possible mass. 

• The most interesting result obtained is that there are possible WDs whose mass 
exceeds the Chandrasekhar limit and is found to be about 2.3M . For instance, for 
tFmax = 20 and B = 2.94 x 10 15 G, we obtain such a WD with p c = 2.2 x 10 9 gm/cc. 
Interestingly, the nature of mass-radius relations does not depend on the value of 
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t-Fmaxi however ep max determines how relativistic the system is. For instance, the 
Chandrasekhar limit is not exceeded for epmax — 2, no matter what p c is. 

• As the magnetic field increases the WDs become more compact in size. 

• The minimum magnetic field required to have a 2.3M WD is B = 2.94 x 10 15 G. 
The magnetic field of the original star of radius R & , which collapses into a WD of 
radius ~ 10 8 cm with the above field at its center, then turns out to be ~ 6 x 10 9 
G, based on the flux freezing theorem. Existence of such stars is not ruled out [14J. 

• In principle, the strong magnetic field causes the pressure to become anisotropic [T5] . 
[16J. As a result, the combined fluid-magnetic medium develops a magnetic tension 
|17j . leading to a deformation in the WDs along the direction of the magnetic field. 
The WDs hence could adopt a flattened shape [15] . 

• The flattening effect due to magnetic field leads to super-Chandrasekhar WDs even 
for smaller magnetic fields. Such relatively weakly magnetized WDs are more prob- 
able in nature. 

• We end by addressing the possible reason for not observing such a high field yet in a 
WD. This could be due to the magnetic screening effects on the surface of the WD 
if it is an accreting one. In this case the current in the accreting plasma depositing 
on the surface, presumably creates an induced magnetic moment of sign opposite 
to that of the original magnetic dipole, thus reducing the surface magnetic field of 
the WD, without affecting the central field. Hence by estimating the surface field 
alone one should not assume the rest. 
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